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SUMMARY 


This paper describes two very important and useful features available in 
Level 17.0 version of NASTRAN. The first one is the APPEND feature which is 
applicable in the case of real eigenvalue analysis. This feature permits the 
addition of new eigenvalues and eigenvectors to those already computed in a 
previously checkpointed run without re-executing the entire problem. The second 
feature is the CONTINUE feature which is applicable in the case of transient 
analysis of coupled equations. This feature enables the integration of coupled 
equations to be continued beyond the last (or from any earlier intermediate) 
output time for which the solution was obtained in a previously checkpointed 
run (without re-executing the entire problem). The paper illustrates the use 
of these two features by suitable examples. 


INTRODUCTION 


The checkpoint/restart feature available in NASTRAN is a very useful capa- 
bility that permits the restarting of previously checkpointed runs without 
re-executing the entire problem. The modules that need to be executed on re- 
start are determined by the nature of modifications made to the checkpointed 
data by the user upon restart. Depending on these modifications, four types 
of restarts may be identified. These are the Unmodified Restart, Modified 
Restart, Rigid Format Switch and Pseudo Modified Restart. The details of these 
restarts are discussed in Reference 1. It is to be noted, however, that the 
user does not explicitly specify the type of restart. It is implied and auto- 
matically determined by the changes made to the NASTRAN data deck upon restart. 

It is important to note that a restart only determines the modules that 
need to be executed during the current run. A restart does not, in general, 
affect the logic of execution within a given module. There are, however, two 
important exceptions to this. These are the APPEND and CONTINUE features 
available in Level 17.0 version of NASTRAN. Thesa are discussed in detail in 
the following sections. 


THE APPEND FEATURE 


In real eigenvalue analysis, it is frequently necessary to add new eigen- 
values and eigenvectors to those already computed in a previous run. The APPEND 
feature available in Level 17.0 version of NASTRAN makes it possible to do this 
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without re-executing the entire problem. It is available only when using the 
Inverse Power, Determinant and Tridiagonal Reduction (FEER) methods of eigenvalue 
extraction in Rigid Formats 3, 10, 11, 12, 13 and 15 (Displacement approach) and 
in Rigids Formats 10 and 11 (Aeroelastic approach). This feature is particular- 
ly valuable in the case of large order eigenvalue problems. 

In-order to use the APPEND feature, the user first requests a checkpoint 
of an eigenvalue problem employing one of the three above-mentioned methods of 
eigenvalue extraction. This run can terminate for any reason so long as the 
READ (Real Eigenvalue Analysis-Displacement approach) module finds at least 
one eigenvalue and one eigenvector and the LAMA (eigenvalue) and PH IA (eigen- 
vector) files are successfully checkpointed. The READ module also sets the 
parameter NEIGV to be equal to the number of eigenvalues and eigenvectors 
found on this checkpoint run. 

The user then restarts and activates the APPEND feature by changing either 
the METH0D card in the Case Control Deck and/or the EIGR card in the Bulk Data 
Deck so as to force the re-execution of the READ module. The method of eigen- 
value extraction used in the restart need not be the same as that used in the 
checkpoint run, but the structural model and the constraint data must be the 
same. Also, the user must ensure that the range of eigenvalues specified on 
the EIGR Bulk Data card for the restart does not include the eigenvalues that 
have previously been found and checkpointed. It is left to the user to satisfy 
this requirement. The program does not check for this condition. 

The APPEND feature causes the READ module to retrieve the eigenvalues and 
eigenvectors from the previously checkpointed LAMA and PHIA files, respectively, 
(this retrieval is done in subroutine READ7 within the READ module; the number 
of eigenvalues and eigenvectors retrieved is indicated by a user information 
message) and to subsequently combine them with the newly computed results. 

This is shown schematically by the flow diagram in Figure 1. The eigenvalues 
and eigenvectors output by the restart include those previously checkpointed. 
Also, the resulting eigenvectors are normalized according to the method of 
normalization specified in the restart. 

In certain cases of restart, the user may not be interested in retrieving 
all the NEIGV eigenvalues and eigenvectors found on a checkpoint run. In such 
cases, the user may retrieve only the first n (n < NEIGV) eigenvalues and eigen- 
vectors from the LAMA and PHIA files, respectively, by resetting the parameter 
NEIGV in the restart to be equal to n (this capability is not available in 
Level 17.0, but will be in Level 17.5.) by means of a PARAM statement just 
before the READ module in the DMAP sequence. This is done by means of a DMAP 
alter in the Executive Control Deck of the restart. 


THE CONTINUE FEATURE 


In transient analysis, it is frequently necessary to continue the integra- 
tion of the coupled equations beyond the last (or from any earlier intermediate) 
output time for which the solution was obtained in a previous run. Thus, the 
initial time for the new run is to be a specified output time of the previous 
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run. Also, the displacements, velocities and accelerations corresponding to the 
specified output time of the previous run are to be used as the initial condi - 
tions for the new run. The CONTINUE feature available in Level 17.0 version of 
NASTRAN makes it possible to do this without re-executing the entire problem. 

It is available In both Rigid Formats 9 and 12 (Displacement approach). This 
feature can be particularly valuable in the case of large order transient analy- 
sis problems involving extended integrations. 

In order to use the CONTINUE feature, the user first requests a checkpoint 
of a coupled transient analysis problem in the normal manner. This run can 
terminate for any reason so long as the TRD (Transient Analysis-Displacement 
approach) module computes the solution for at least one output time and the 
UDVT (displacement-velocity-acceleration) file is successfully checkpointed. 

The T0L (list of output times) file would have been previously checkpointed 
subsequent to the execution of the TRLG (Transient Load Generator) module. 

The TRD module also sets the parameter NC0L to be equal to the number of output 
time steps (which is also equal to one-third the number of columns in the UDVT 
matrix) in the checkpoint run. 

The user then restarts and activates the CONTINUE feature by changing any 
one or more of several cards either in the Case Control Deck (DL0AD, N0NLINEAR, 
TSTEP cards) and/or in the Bulk Data Deck (TSTEP, DAREA, DL0AD, F0RCE, etc.) 
that define either the dynamic loading and/or the time step selection. This 
forces the re-execution of both the TRLG and TRD modules. The dynamic loading 
and/or the time step selection in the restart need not be the same as that used 
in the checkpoint run, but the structural model and the constraint data must 
be the same. 

The CONTINUE feature causes the integration of the coupled equations to 
continue from a specific initial time. For normal restarts (in which the 
checkpointed value of the parameter NC0L is automatically used), this initial 
time is the last output time of the checkpoint run. However, in certain cases, 
the user may wish to restart the integration from an intermediate (rather than 
from the last) output time of the checkpoint run. In such cases, the user 
should reset the parameter NC0L to correspond to the desired intermediate out- 
put time by means of a PARAM statement just before the TRLG module in the DMAP 
sequence. This is done by means of a DMAP alter in the Executive Control Deck 
of the restart. 

The output of the restart does not include the solutions of the checkpoint 
run, but only those solutions that are computed by the restart. Also, any ini- 
tial conditions specified in the restart data are ignored since the solution is 
continued by using the displacements, velocities and accelerations corresponding 
to the specified output time of the checkpoint run as initial conditions for 
the restart. 

The first displacement {u^} of the restart (or CONTINUEd run) is given by 
[D] {Ui } = j {?_■] + Pg + P^} + {Ng> + [C] {u n > + [E] , (1) 
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where the matrices [0], [C] and [E] are given by 


CD] » (— [HI + -L- [B] + i [K].) , (2) 

\At 2 2 At J / 

[C] " [Ml - 1 [K]j , (3) 

and [El = [M] + -J- [B] - i [K]l (4) 

At 2 2 At J 

and {P 0 } = [K] {u n } + [B] {u p } + [M] {u p } , (5) 

..2 

(u_-j } - {u n J - At {G n } + oy , (6) 

{u^> = <u n } - {U n > At , (7) 

{P_.j} = [M] {u‘ n } + [B] {u_-j } + [K] (u_ 1 } . (8) 


In the above equations, [Ml, [B] and [Kl are the mass, damping and stiffness 
matrices, respectively; {u n >, (%} and {ti n } are the displacements, velocities 
and accelerations, respectively, at the specified output time t p of the check- 
point run, and At is the initial time step for the restart. {P^} is the load at 
time t = t n + At and {Nq > is the initial non-linear load. 

The assumptions represented by Equations (6) through (8) introduce errors 
in the restart. These inherent errors may be minimized by selecting the initial 
time step in the restart to be the same as the time step used in the checkpoint 
run gust before the restart. 


EXAMPLES 


In order to illustrate the use of the APPEND and CONTINUE features dis- 
cussed above, two examples were selected. Example 1 is the analysis of the 
transverse vibrations of a 10-cell beam. The finite element model used is 
shown in Figure 2. Example 2 is the transient analysis of a 1000-cell string 
(travelling wave problem). This is the same as NASTRAN Demonstration Problem 
No. 9-2-1 (Reference 2). The finite element model used is shown in Figure 3. 
Both these problems were run on the CDC CYBER computer using a post-Level 17.0 
version of NASTRAN. 
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Discussion of Example 1 


All the twenty (20) natural frequencies of the model of Example 1 were 
first computed. These are listed in Table 1 to facilitate comparison with 
subsequent results. In order to illustrate the use of the APPEND feature, a 
checkpoint run was then made using the Inverse Power method and the first eight 
modes were computed, These are presented in Table 2. Using the results of 
this checkpoint run, restarts were then made under various conditions using dif- 
ferent extraction methods. The results obtained are presented in Table 3. As 
can be seen, all these restarts involve the retrieval of the eight modes of the 
checkpoint run of Table 2. 


Discussion of Example 2 


A transient analysis of the model of Example 2 was first made using twenty 
(20) time steps. (All runs for this example were made on Rigid Format 9 using 
time steps of 0.0005 seconds.) The displacements for point 10 are listed in 
Table 4 to facilitate comparison with subsequent results. In order to illus- 
trate the CONTINUE feature, a checkpoint rur, was then made using only ten (10) 
time steps. The displacements for point 10 obtained in this run are presented 
in Table 5. As can be seen, these results are merely a subset of those in Table 
4. The integration of Table 5 was then CONTINUEd for ten (10) more time steps 
by restarting from the last output time in Table 5 (0.005 second). The results 
are presented in Table 6. In order to illustrate the resetting of the parameter 
NC0L, an additional restart was made by changing the value of NC0L to 6 (from 
its original value of 11) and the integration of Table 5 was CONTINUEd for 
fifteen (15) more time steps. This thus involved starting from an initial time 
of 0.0025 second. The results of this run are shown in Table 7. A comparison 
of the results of Tables 6 and 7 with the corresponding results in Table 4 
reveals the inherent. errors caused by the CONTINUE feature. 


SUMMARY AND CONCLUSIONS 


Two very important and useful features available in Level 17.0 version of 
NASTRAN have been described. The first one is the APPEND feature which is ap- 
plicable in the case of real eigenvalue analysis. This feature permits the 
addition of new eigenvalues and eigenvectors to those already computed in a 
previously checkpointed run without re-executing the entire problem. The 
second feature is the CONTINUE feature which is applicable in the case of tran- 
sient analysis of coupled equations. This enables the integration of coupled 
equations to be continued beyond the last (or from any earlier intermediate) 
output time for which the solution was obtained in a previously checkpointed 
run (without re-executing the entire problem). The use of these two features 
has been illustrated by means of two suitable examples. 
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Table 1. Results for Example 1 
(using Givens method) 


Mode 

no. 

Natural frequency 

(Hz) 

l’ 

1.591560E-02 

2 

6.366B79E-02 

3 

1.433160E-01 

4 

2.550698E-01 

5 

3.994578E-01 

G 

5.7750BOE-01 

7 

7. 9091 6 IE-01 

G 

1.042057E+00 

9 

1 ,33276'lE-i'OO 

10 

1 . 76G489E+00 

11 

2.091G90E+00 

12 

2.549921 [+00 

13 

3.0B6530E+00 

14 

3.71 2652E+00 

15 

4.440216E+00 

16 

5.273030E+00 

17 

6.187597E+00 

IS 

7.09701 3E+00 

19 

7.814544E+00 

20 

8.095071 E+OG 


Table. 2. Results of Checkpoint Run of Example 1 
(using Inverse Power method) 


Mode no. 

Natural frequency 

in Table 1 

(Hz) 

1 

1.S91560E-02 

2 

6.366879E-02 

3 

1 . 433160E-01 

4 

2. 550698E-01 

5 

3.99 4578E-01 

6 

5.775080E-01 

7 

7.909161E-01 

8 

1.042057E+00 

Eigenvalue extraction data 

FI = 0.0 Hz 

NE - 8 

F2 = 1.2 Hz 

ND = 8 
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Table 3. Results of Restart Runs of Example 1 Using Checkpoint Run of Table 2 


Invers 

e Power method 

Determinant method 

Tridiagonal Reduction 
{ PEER) method 

Remarks 

Mode no. 
in Table 1 

flatural frequency 
(Hz) 

Mode no. 
in Table 1 

Hatural frequency 
(Hz) 

Mode no. 
in Table 1 

Hatural frequency 
(Hz) 

1 

1.591560E-02 

1 

1 .59156OE-02 

1 

1 . 591 560E-02 

* 

2 

6.366379E-0Z 

2 

6.356379E-02 

2 

6.366B79E-02 


3 

1 .4331 60E-01 

3 

1.433160E-01 

3 

1.433160E-01 

Results retrieved 

4 

2.550698E-01 

4 

2.550593E-01 

4 

Z.550693E-01 

from the 

5 

3.994578E-01 

5 

3.99457SE-01 

5 

3.994578E-01 

checkpoint run 

6 

5. 775080E-01 

6 

5.775030E-01 

6 

5.775080E-01 

of Table 2 

7 

7.9091 61 E-01 

7 

7. 9091 61 E-01 

7 

7.9091 "IE-01 


3 

1.0420S7E+00 

8 

1.042057E+00 

3 

1.042057E+00 


12 

2. 549921 E+OO 

10 

1.766489E-f-00 

g 

1.332764C+00 


13 

3.086530E+00 

11 

2.091890E+00 

10 

1.766489E+00 


14 

3.71 2652E+0Q 

12 

2. 549921 E+00 

n 

2.091390E+00 


15 

4.44O2I6E+O0 

13 

3.085530E+00 

12 

2. 549921 E+00 


16 

5.273030E+00 

14 

3. 7 12652 E+00 

13 

3 . 086530E+00 


17 

6.1B7597E+00 

15 

4.44021 6E+00 

14 

3.71 2652E+00 

Results computed 

18 

7.09701 3E+00 



15 

4.44021 6E+00 

by the 

19 

7.314544E+00 



16 

5 . 273030E+00 

method selected 





17 

6.187597E+00 






18 

7.09701 3E+00 






19 

7.614544E+00 






20 

3.095071 E+00 


Eigenvalue extraction data 


F! = 3.0 Hz 

00 

II 

UJ 

FI = 2.0 Hz 

HE = 6 

FI = 5.0 Hz 

HD = 11 


F2 = 9.0 Hz 

HD = 16 

F2 = 6.0 Hz 

HO = 14 






Table <1. Results for Example 2 


°kk;w A |, 
i00, i QtMfJT 


IS 

Y 


Time 
(sec. ) 

— 

Time step 
no. 

Displacement 
of f’oint 10 

0.0 

1 

1 , OOOOOOF.tOO 

0,0005 

2 

1 .797927 0 00 

0.0010 

a 

1.7071022*00 

0.0015 


1 .755303!' rOQ 

0.0020 

5 

1.0906052+00 

0.0025 

6 

1.5099O7C+OQ 

0.0030 

7 

1 .-164-1972+00 

0.0035 

fl 

1. 334031IM0Q 

0.0040 

9 

1. 2203-IRE 100 

0. 00*10 

10 

1 .1279215+00 

0.0050 

11 

1 , 0-196391: tOO 

0.0055 

12 

9. 691 G4 52-01 

0.0060 

13 

0.733-1-142-01 

0.0060 

14 

7.607323E-Q1 

0.0070 

15 

G. 4161592-01 

0.0075 

16 

5. 3Q04G0E-01 

0.0000 

17 

-1. 330371 E-Ql 

0.0005 

ia 

3.4 95874 C— 01 

0.0090 

19 

R. 6595-11 E-01 

0.0095 

20 

1. 72-14 13E-01 

0.0100 

21 

6,72gmm:-02 


Table 5. Results of Checkpoint Rim of Example R 


Time 
(sec. ) 

~ — ■ ""1 

Time step no, 
of Table -1 

Displacement 
of Paint 10 

0.0 

1 

i .aooooon+oo 

0.0005 

2 

1. 797927 E+QO 

0,0010 

3 

1.7071OHE+OO 

0.0015 

-1 

1.7553032+00 

0.0020 

5 

1.6906052+00 

0.00?5 

6 

1.5899872+00 

0.0030 

7 

1 ,4 6-1 -1971; tOO 

0.0035 

0 

1,3348812+00 

0.004O 

9 

1 . 2203-122+00 

0.00 

10 

1.1279212+00 

O.OOi 

11 

1 .0-196392400 
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Table 6. Results of Restart Run of Example 2 
Using Checkpoint Run of Table !> 


Time 

(sec.) 

... 

Time step no. 
of Table 4 

Displacement 
of Point 10 

0.0050 

11 

1.049639E+00 

0.0055 

12 

9.644270E-01 

0.0060 

13 

8.67398GE-01 

0.0065 • 

14 

7.575025E-01 

C.0070 

15 

6.432689E-01 

0.0075 

16 

5.357538E-01 

0.0080 

17 

4.405755E-01 

0.0085 

18 

3, 543471 E-01 

0.0090 

19 

2.677839E-01 

0.0095 

20 

1. 72941 3E-01 

0.0100 

21 

6.924284E-02 


Table 7. Results of Restart Run of Example 2 Using Checkpoint Run 
of Table 5 and With Parameter NC0L Reset to 6 


Time 

(sec.) 

Time step no. 
of Table 4 

Displacement 
of Point 10 

0.0025 

6 

1.589987E+00 

0.0030 

7 

1.470335E+00 

0.0035 

8 

1 .344651 E-t-OQ 

0.0040 

9 

1.229273E+00 

0.0045 

10 

1.132560E+00 

0.0050 

11 

1.049933E+00 

0.0055 

' 12 

9.681 384E-01 

"0.0060 

13 

8.746610E-01 

0.0065 

•14 

7.659208E-01 

0.0070 

15 

6.489350E-01 

O.OG75 

16 

5.357538E-01 

0.0080 

17 

4.349093E-01 

0.0085 

18 

3.459289E-01 

0.0090 

19 

2.605215E-01 

0.0095 

20 

1.692299E-01 

0,0100 

21 

6.359591E-02 
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Figure 1. Flow Diagram for the Real Eigenvalue Analysis Module, READ. 
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Parameters : 

Length of each BAR element = l = 3 1000.0 In. (25.4 m) 

Young's modulus 3 Shear modulus * 1.0 x 10 7 psi (6.89477 x 10 10 fj/m Z ) 

Hass density = 0.9869604 lb.-sec 2 /fn. 4 (1.05475 x 10 7 Kg/m 3 ) 

Area of cross section of each BAR element = 9,869604 in. 2 

3 6.36747 x 70" 3 m Z 

Area moments of inertia 3 1.0 x 10 6 in. 4 

3 0.41623 m 4 


Figure 2. Representation of 10-Cell Beam of Example 1 
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Finite Element Model 


Parameters : 

= ^- = 10 7 (stiffness units) for all i 
= yAx = 10 (mass units) for all 1 
Loading : 

The initial displacements are given by: 
u i = 0.2 (1-1) for 2 i i ^.11 

u j = 0.2 (21-i ) for 11 < i < 21 
and 

= 0.0 for i >.21 

Figure 3. Representation of 1000-Cell String of Example 2 
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